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ABSTRACT 


This  research  effort  investigates  system  requirements,  configu¬ 
rations,  and  associated  performance  characteristics  of  maximum- 
likelihood  space-diversity  receiving  systems.  From  the  maximum 
likelihood  principle,  and  the  multidimensional  Karhunen-Loeve  expan¬ 
sion,  the  continuous  decision  equations  for  space-time  detection  of 
noisy,  random  multipath  signals  are  developed.  Interpretation  of 
these  equations  yields  the  optimum  space-time  receiving  system  con¬ 
figuration  requirements. 

The  performance  characteristics  of  the  resultant  systems  are  then 
analyzed  ’o  determine  tfce  effects  of  channel  multipath  structure 
(multipath  delay  and  power  division  among  the  paths),  space-time 
correlation  properties  of  the  incident  processes,  and  the  specific  tem¬ 
poral  correlation  introduced  by  the  array  geometry.  It  is  shown  by  a 
series  of  case  studies,  that  for  both  the  single  element  coupling  as  well  as 
array  multiple  element  coupling,  that  increasing  the  multipath  delay  factor 
results  in  decreased  system  performance  capability  for  fixed  power  of 
the  signal  and  noise  processes.  Similarly,  the  performance  capacity  is 
degraded  as  the  available  signal  power  tends  to  distribute  more  equally 
over  the  multiple  transmission  paths.  These  effects  are  attributed  to  the 
loss  of  effective  signal  energy  concentration,  resulting  in  a  lower  effective 
detectability  signal-to-noise  ratio.  An  investigation  of  these  effects  upon 
system  performance,  due  to  array  geometry  (receiving  element  spacing) 

shows  that  performance  is  enhanced  by  increasing  the  separation  distance 

ii 


of  the  elements  for  a  given  multipath  situation.  This  phenomena  is 
the  result  of  a  more  compatible  array  beam  pattern  that  comes  about 
from  increased  spacing.  In.  effect,  with  increased  spacing,  the  main 
lobe  of  the  pattern  is  narrowed,  while  the  side  lobes  are  optimally 
suppressed  by  the  optimal  cross- coupling  of  elements.  Finally  the 
analysis  illustrates  how  optimum  space-time  receiving  systems 
capability  results  from  a  joint  consideration  of  the  coupling  and  tem¬ 
poral  processing,  rather  than  aidichdtoniiaat'idn-of -thelprelildmainto  a 
spatial  (antenna)  problem  and  a  (temporal)  data  processing  problem. 

As  such,  the  analysis  yields  a  more  concise,  broader,  interpretation 
of  system  design  requirements  and  associated  performance  character¬ 


istics. 
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CHAPTER  I 


INTRODUCTION 

In  most  modern  communications/radar/sonar/seismology 

receiving  systems,  the  concept  of  space  diversity,  or  multi-element 

12  3 

array  reception,  is  of  considerable  importance.  *  ’  In  such  a  system, 
the  receiver  consists  of  an  array  of  receiving  elements,  antennas,  hydro 
phones,  or  gcophones  distributed  in  space.  These  elements  couple  the 
data  processor  to  the  electromagnetic  or  acoustic  medium  containing  the 
desired  information  processes.  It  is  through  these  spatially  distributed 
elements  that  the  spatio-temporal  (space-time)  wave  function  processes 
are  transformed  into  temporal  (time)  voltage  waveforms.  Following  the 
array  elements  is  the  receiver  stage,  which  steers  the  beam  pattern  of 
the  array  in  the  direction  of  the  incident  information  processes  thus 
performing  "beam  steering"  and  "phasing"  of  the  array.  This  stage, 
in  addition  to  the  array  geometry,  attempts  to  optimally  exploit  the 
spatial  structure  of  the  incident  process  (analogous  to  the  temporal 
processing  in  the  frequency  filtering  and  combining  stages  of  the  data 
processor).  Thus  the  array  and  its  associated  steering  stages  can  be 
considered  to  be  a  spatial  filter  as  well  as  a  temporal  data  processor. 

Following  the  spatial  filter  is  the  temporal  data  processor;  it 
is  here  that  the  receiver  performs  operations  (prediction,  smoothing, 
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correlation,  etc. )  on  the  spatially  filtered  data  to  carry  out  the  desired 
detection,  classification,  and  estimation  functions  that  yield  the 
desired  information  (signal,  no- signal;  target,  no-target;  range,  doppler, 
etc. ). 


Thus  the  total  receiver  structure  shown  in  Fig.  1-1  consists  of 

a  spatial  filter  and  a  temporal  filter.  In  order  to  attain  ''optimum" 

reception,  one  must  treat  the  spatio-temporal  processing  jointly.  Previous 

work  has  partitioned  the  over-all  reception  problem  into  a  spatial  solution 

4  5  6 

(antenna  array  design,  beam  pattern  synthesis)  *  *  and  a  temporal 

solution  (classical  statistical  communication  theory,  detection  theory, 

etc.).  Recent  studies  have  shown  the 

utility  of  working  with  the  entire  receiver  system  as  a  unit  by  considering 

the  combined  space-time  receiving  system  operations  to  obtain  optimum 

receiving  systems  with  improv\  J  performance  characteristics^ 

The  temporal  data  processor  can  be  either  an  analog  system 

operating  upon  continuous  waveforms  or  a  discrete  system  in  which  the 

coupled  processes  are  sampled  (A/D  conversion)  at  some  stage  and  pro- 

25 

cessed  with  a  digital  filter  (digital  computer).  The  analysis  that  follows 
treats  the  waveform  case  and  develops  the  expressions  for  the  optimum 
continuous  decision  equations  or  test  statistic  computers.  Interpretation 
of  these  equations  yields  the  resulting  optimum  continuous,  analog 
system  configuration.  This  structure  is  highly  dependent  upon  the  second- 
order  spatio-temporal  properties  of  the  incident  processes,  the  spatial 


19,  20,  21,  22, 
23,  24 
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filtering  introduced  by  the  array  geometry,  and  the  temporal  filtering 
of  the'frequency  filters,  combining  sections  and  beam  forming  stages. 

In  this  research  a  model  will  first  be  developed  for  the  informa¬ 
tion  channel  physics  to  which  the  receiver  will  be  coupled.  From  this 

model  the  problem  structure  is  formulated.  Then,  using  the  maximum 

11  12  26  27 

likelihood  principle,  *  '  1  the  coupled  process  expansion  in.  orthog¬ 

onal  coordinates  or  sufficient  statistics,  and  the  Karhunen-Loeve  series 
expansion  for  random  processes,  the  space-time  decision  equations  are 
developed.  The  performance  characteristics  of  the  resulting  systems 
(detectability,  and  receiver  operation  characteristics)  are  then 
developed,  analyzed,  and  interpreted.  Finally,  the  performance  for  an 
array  receiving  system  in  a  random  multipath  environment  is  presented 
to  illustrate  the  effects  of  array  geometry,  the  signal  multipath  structure, 
and  the  spatio-temporal  properties  of  the  incident  random  processes  upon 
system  performance.  From  such  an  analysis  one  obtains  guidelines 
which  enhance  the  design  and  implementation  of  more  optimum  detection 
estimation,  classification  and  discrimination  systems.  The  analysis  and 
results  are  found  in  terms  of  process  intensity,  power,  power  spectra, 
array  geometry,  etc.  Thus,  this  paper  proposes  to  resolve,  in  terms 
of  the  parameters  of  physics  (time,  energy,  and  geometry)  the  problems: 

1.  For  a  given  space-time  structure  of  signal  and  inter¬ 
ference  processes,  what  is  the  optimum  system  configuration 
in  terms  of  geometry  and  temporal  structure? 


For  such  systems,  what  are  the  expected  optimum  and 
actual  (mis-matched)  system  performance  character- 
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SYSTEM  ENVIRONMENT  AND  PROCESS  MODEL 


The  waveform  processes  incident  upon  the  array  of  receiving 
elements  consist  of  both  noise  and  signal  components.  The  noise 
components  are  due  to  a  combination  of  all  random  interfering  processes 
that  exist  within. the  receiving  environment.  Typical  examples  in  the 
sonar  case  include  sea  noise,  noises  due  to  fish,  noises  due  to  other 
ships  and  submarines.  The  signal  components  are  due  to  signals  trans¬ 
mitted  by  the  source,  as  in  point-to-point  communications,  signals 
reflected  by  the  target,  as  in  active  radar,  sonar,  seismology,  or 
random  signals  propagated  by  the  source  as  in  passive  radar /sonar  or 
radio  astronomy. 

These  signals  can  be  either  random,  deterministic,  or  either 

with  random  or  unknown  parameters?*'  This  paper  will  treat 

the  purely  random  signal  case  in  which  the  incident  signal  components 

are  assumed  to  be  sample  functions  from  a  zero-mean  (no  specular 

component),  gaussian  process,  produced  by  a  channel  or  cross  section 

containing  multipath  and  a  scattering  mechanism.  Such  a  random, 

multipath,  information  signal  is  commonly  encountered  in  tropo-scatter 

communication  links,  ^  telemetry  links  over  randomly  fading 

channels?^’  and  in  sonar/radar  detection-estimation  media?’  ^ 

The  multipath  structure  can  arise  from  transmission  through  layered 
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media,  as  in  the  atmosphere  in  radar,  or  underground  earth- rock 

layers  as  in  seismology,  in  which  parts  of  the  signal  are  reflected  and 

29 

take  different  paths  to  the  receiver  array.  The  random  fading  can 

occur  in  channels  in  which  the  transmission  media  is  moving  or  varying 

in  density,  such  as  transmission  of  acoustic  signals  through  the  sea, 

28 

etc.  Thus  the  received  signal  waveform  can  be  classified  as  a  random 

signal  or  a  randomly  faded,  deterministic  signal  that  has  a  scattering, 

21,  22,  24 

multipath  structure  due  to  the  channel  media.  As  such,  this 

paper  deals  with  the  space-time  (array)  detection  of  a  random,  zero- 
mean,  multipath  signal  in  additive  gaussian  noise.  The  signal  will  be 
taken  as  purely  random  as  is  often  encountered  in  passive  radar,  sonar, 
and  radio  astronomy. 

2-1.  Multipath  Structure 

The  multipath  structure  will  be  assumed  resolvable  and  charac- 
178  30  31 

terized  by  delay  A. '  *  *  ’  This  delay  is  due  to  the  difference  in 

transmission  path  length  between  source  and  receiver;  it  is  assumed 
that  the  signals  of  each  path  arrive  from  the  same  direction  with  devia¬ 
tions  from  the  angle  of  arrival  small.  The  signal  can  be  expressed  as 


Zch(t)  =  I  zchm<t-mA) 


(2.1) 


where  p  is  the  number  of  paths  and  A  is  the  characteristic  delay. 
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i 

2-2.  Channel  Process  Model 

i 

The  channel  noise  process  is  assumed  taken  as  a  zero-mean, 
gaussian  process  that  additively  masks  the  random  signal  waveform. 

The  total  incident  process  is 

ria(t>  ■  2ch(t)  +  <2'2> 

The  random  processes  described  above  are  incident  upon  an 
array  of  K  elements  distributed  in  space.  The  geometry  is  shown  in 
Fig.  2-1.  The  location  of  each  element  relative  to  a  known  co-ordinate 
system  is  described  by  its  respective  position  vector  j>.  The  ith  element 
is  located  in  space  at  a  point  described  by  the  vector  p.. 

T 

•  Pi  =  CPxiPyiPki^  (2-3) 

The  incident  wave  interacts  with  the  receiver  elements  to  produce  at  the 
output  terminals  of  each  element  a  received  process  r(t).  The  output 
of  the  ith  element  is  r.(t).  For  a  K  element  array,  the  set  of  output 
processes  are  written  as  the  vector  process 

r(t)  =  [^(t)  r2(t)  .  . .  rk(t)]T  (2.4) 

It  is  assumed  that  the  process  r^(t)  results  from  passage  of  the  channel 
process  through  the  receiver  element  with  linear  impulse  response 

i 

h£(t,  u),  allowing  for  imperfect  receiving  elements.  The  received  process 


is  thus 


Figure  2-1..  Arr^ay  Geometry 
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t 

r.(t)  =  hjt,  u)  rch(u)du  (2.5) 

Each  received  waveform  consists  of  a  signal  process  plus  a 
noise  process.  In  addition  there  may  be  noise  due  to  the  receiver  element 
or  receiver  "front  end,"  which  can  be  included  as  an  additional  post 
additive,  coupling  system  noise  process  n  (t).  Thus,  the  total  received 
process  at  the  output  of  the  ith  element  is 

t 

ri(t)  =  ^  he(t,  u)  r^tu)  du  +  nr.(t) 
t 

o 

=  z.(t)  +  n.(t)  (2.  6) 

where  z.(t)  is  the  received  signal  component  and  the  two  noise  processes 
can  be  added  to  produce  the  single  noise  interference  n.(t).  The  set  of 
received  processes  given  in  (2.  4)  is  now 

r(t)  =  z(t)  4  n(t)  (2.7) 


For  the  signal  and  noise  independent,  zero-mean,  gaussian 
processes,  the  covariance  function  matrix  of  r(t)  is  defined  as 


with 


Kr(t,  u)  = 

E[r(t)  r"(u)] 

(2.8) 

Kr(t,  u)  = 

K  (t,  u)  +  K  (t,u) 

^  n 

(2.9) 

where 
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K  (t,  u)  =  E[z(t)  z  (u) ]  (2.10) 

2  1  ' 

and 

Kh(t,  u)  =  E[  n(t)  n  (u)] 

N 

=  6(t-u)  I  +  K  (t,  u)  (2.11) 

L  —  c 

The  elements  of  the  process  covariance  matrices  can  be  termed 
self- correlation  functions  for  elements  on  the  diagonal  of  the  covariance 
matrix  and  the  off  diagonal  terms  are  a  measure  of  the  correlation  in  the 
processes  between  elements.  Thus  the  correlation  functions  K  (t,  u) 

ij 

can  be  denoted  "self-link"  correlation  functions  if  i  -  j  and  "cross  link" 

functions  if  i  i  j.  The  cross-link  correlation  functions  are  directly 

related  to  the  spatial  filtering  produced  by  the  array  geometry  in  terms 

of  the  spatial  correlation  structure  of  r(t).  This  dependence  results 

from  the  functions  being  dependent  on  the  temporal  variables  and  spatial 

variables  (spatio-temporal  correlation  properties)  of  the  incident  process 

and  upon  the  spatial  variables  of  the  array  geometry. 

The  dependence  upon  array  geometry  can  be  seen  by  examining 

the  received  signal  waveform  vector  r(t).  Since  the  signal  is  assumed 

to  be  a  plane  wave  with  incidence  angle  unit  vector  pg ,  the  wavefront 

6  8 

is  time  delayed  due  to  spatial  distribution  of  the  receiver  elements.  ’ 

The  received  process  in  the  ith  element,  located  by  the  position  vector 
p.  is  delayed  by  the  amount  t\  .  The  received  signal  vector  is  written 


z(t)  = 


z2(t) 


zk(t) 


2  (t  -  Tx) 


2  <t-T2) 


z  (t  -  Tk) 


(2.12) 


where  t.  is  given  by  the  known  wave  front  delay  equation 


4,  5,6 


T.  = 
1 


P  ’  P- 

rS  1 


(2.13) 


with  v  the  medium  propagation  velocity. 
Using  (2.12)  and  (2.10)  gives 


K  (t,  u)  =  [Kz  (t- t.s  u- t.)1 

Z  XJ  1  J 


for  all  i,  j  £  K 


(2.14) 


The  same  analysis  is  valid  for  the  noise  correlation  matrix.  Thus 

the  covariance  matrices  are  functions  of  temporal  and  spatial  variables 

and  as  such  are  proportional  to  the  original  spatio-temporal  correlation 

functions.  The  general  structure  of  such  functions  in  terms  of  medium 

18  32 

and  aperture  geometry  have  recently  been  investigated  by  Middleton.  * 
The  effects  of  the  wave  front  delay  of  the  signal  components  is 
counteracted  by  "steering"  the  array  beam  pattern.  The  array  steering 
process  is  accomplished  by  passing  each  received  process  through  an 
appropriate  delay  given  by  the  negative  of  the  wave  front  delay.  Thus,  the 
ith  process,  after  steering,  becomes,  the  prime  denoting  the  steered 


process 
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;  r f 

r!(t)  =  V  r.(u)  6(u-t-T-)du  =  r.(t-M\) 


=  Z.{t-T.-rT.)  +  n.(t+T.) 

1  1  1  l  1 


(2.15) 


=  z(t)  -5-  n.(t+T.) 

l  l 


z  (t)  -r  n.T(t) 
i  i 


Similarly  the  jth  steered  process  is 


rf'(t)  =  z.(t)  -f  n.(t-fT.)  =  z(t)  -i-  n?(t) 
3  J  J  J  1 


(2.16) 


The  steering  vector  process  is 


L^t).  =  ^  Hs(t,u)r(u)du 


where  the  steering  matrix  Hs(t,  u)  is  a  linear  system  with  matrix 


impulse  response 


H  (t,  u)  =  [6(u-t-r.)  6..  ] 
^  1  1J 


The  steered  signal  vector  becomes,  for  homogeneous  signal  process 


z‘(t)  =  z(t)  1 


(2.17) 


n'(t)  =  [^(t+Tj)  n2(t+T2)-->  ^(t-f Tk) ] 

The  steering  section  of  the  receiver  aligns  the  signal  components 
of  the  received  process  and  introduces  a  time  shift  in  the  received 


noise  processes. 


structures  will  reflect  these  new  properties  and  show  the  effects  of 
altering  the  noise  field  space-time  properties. 
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2-3.  Hypothesis  Structure 

The  hypothesis  structure  to  be  considered  is  the  binary  signal; 
no-signal  case.  Under  hypothesis  the  received  process  is  the  sum 
of  signal  and  noise  components 

r(t)  =  z(t)  +  n(t) 

under  H  the  received  process  is  due  to  noise  alone,  thus 
o 

r(t)  =  n(t) 

The  correlation  function  matrix  of  r(t)  under  H  is  denoted 

o 

K  (t,  u)  =  K  (t,  u)  (2.19) 

£  n 

and  under  H^,  the  correlation  matrix  is 

K.(t,u)  =  K  (t,u)  +  K  (t,u)  (2.20) 

1  z  n 

The  above  model  and  hypothesis  structure  define  the  physical  channel 
that  is  to  be  studied.  The  channel  model  is  general  enough  to  reflect 
the  physical  properties  of  a  wide  range  of  important  information 
channels  of  technological  interest.  ’  *  The  analysis  of  the  detection 

problem  receiver  structure  s  and  system  performance  characteristics 


will  reflect  the  model  and  hypothesis  structure  developed  above.  As 
each,  the  purpose  of  the  research  effort  defined  at  the  close  of  Chapter  I 
will  be  realized  as  the  optimum  array  (space-time)  diversity  detection 


I 


of  noisy  random,  zero-mean,  multipath,  gaussian  signals.  The  follow¬ 
ing  chapter  will  derive  the  optimum  decision  equation  for  the  channel 


model  above. 


CHAPTER  III 


OPTIMUM  DECISION  EQUATIONS  AND 
SYSTEM  CONFIGURATION 


The  following  chapter,  using  the  previous  channel  model,  estab¬ 
lishes  the  optimum  continuous  likelihood  ratio  test  (LRT/j-decision 
equation  for  the  array/diversity  detection  of  random  signals.  After  this 
relationship  is  established,  the  first  and  second  moments  (mean  and 
variance)  of  the  test  statistic  are  found  for  stationary  and  non- stationary 
random  processes.  Interpretation  of  the  LRT/decision  equation  yields 
the  optimal  array  receiver  structure.  Several  structures  will  be 
investigated. 

3- 1 .  Derivation  of  Decision  Equations 

Using  the  maximum  likelihood  principle,^'  ^  ^  the  expansion 
in  independent  coordinates  of  the  likelihood  ratio  function  defined  as 
is  initiated. 


pMt)|H) 

■  ifegmu  <  \:  Hi 


(3.1) 


For  analysis  purposes,  let  each  of  the  above  densities  be  the  ratio  of 
two  densities,  one  of  the  density  under  Hj  due  non-white  noise,  and 
the  other  the  density  for  r(t)  white  noise  alone.  Thus  (3,1)  becomes 
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L 


The  denominator  is 


r.  2 


p(e I  »0)  iSi „T  expC'2Ci0 )  ) 

,  V  Oj  1 


PfelHJ  ji  «p(-t(A  )2) 

v2tt  ax  awj 


(3.6) 


where 


N  .  N 

.  Z  O.C  1  o 

CT rri  -  X.  7  .  T-  X.  —  X.  4*  . 

■*X.  i  2  i  1  2 

N  N 

o  c  o  o 
Onr  =  ~ r-  -f  X.  =  X.  t  ~=— 

Ao-  2  i  i  2 


(3.7) 


ctt  ~ 

*  wj^  2 

2  Q 

and  X.  and  X.  are  the  eigenvalues  of  the  signal  and  noise  processes 
respectively.  Truncating  the  infinite  product  at  n  and  taking  natural 
logarithms  gives 


ot  _  2 

Aw;  n  r, 

*•  1  tn  s  •! 


L  =  Ln  A 
n 


(rn>  *  {£  Ln  ^7  -  il  (oiT2  -  ^-2)} 

i=l  h  i=l  h  wi 


(3.8) 


■Wj  1=1  J-oi 


Using  the  vector  Karhunen-Loeve  expansion  given  in  (3.4),  and  the 

eigenvalues  defined  in  (3.  7),  the  first-bracket  term  (L^  )  in  (3.  8)  is 

n 


n  ^1  ^ 

i  rr 


L1„  =  I  - 1  $  *<*>  I  — N - —  ds) 

1=1  lT  tT  S  T,  i=l  +  X.1 


(u)du  dt 


Tf  h  tjW  ?iT(“) 

*s*»r 

Ti 


N 

i=l  _2 
2 


r(u)  dudt 


(3.9) 


The  second  term  is  found  in  a  similar  manner  to  be 

'n 


2r 

^  +  X?  ,  Tf  n  ti»  tfa 
2  i  -\  lrr  T, 


^ ■  i  k-v1) n  1  - ^>dtdu 
1=1  —  T.  1=1  +  \ 

i  2  i 


*S&I- 

t  i=l 


n  ^(0  *?(«) 


N 


r(u)  dtdu 


(3.10) 


Ti 


Using  (3.9)  and  (3.10),  taking  the  limit  of  (3.  8)  as  n,  the  number  of  co¬ 
ordinates  approaches  infinity,  combining  terms,  and  putting  constant 
terms  on  the  right  side  of  the  equation  gives 


L  = 


Lim  Ln 
n-*-  <» 


S  rT^[Ml^’  u)  ~  6(t-u)I  jr(u)  dtdu 

T.  ”  ° 


+  \  J  [Mo(t’  U)  "  W~  6(t>u)  :]^  dt  du  J  J  |  h1q 
T.  ”  ° 


(3.11) 
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C3  03 

A  o  2  1  p  2  o 

where  n'  =  L._r>  r  >  (1-i-rr'  X. )  -  )  (1  t  —  X.  )  is  the  new  threshold 
1  n  L,  N  i  N  i 


and  M^(t,  u)  and  M^ft,  u)  are  solutions  to  the  integral  equations 


x. 

^  Kj(t,  u)  M^(u,  z)du  =  6(t-z)£ 


C  K  (t,  u)  M  (u,  z)du  =  6(t-z)I 

J  £  £  — 


(3.12) 


and  defined 


,  *-(t)  ,T(u) 

¥•^14=- 

i=3  -°+XX 

2* 

i 

Mo(t>  U)  =  l  ~W 

i=l  +  x1 

2  i 


(3.13) 


7  8  9 

Following  the  method  of  Price  ’  ’  for  the  single  sensor  case,  by- 


defining 


Mr(t,  u)  =  [6(t-u)I  -  Wx(t,  u)  ] 


3(t,u)  =  ^-[6{t-u)I-  WQ(t,u)] 


(3.14) 


it  can  be  shown  that  the  matrix  filters  W  (t,  u),  W  (t,  u)  are  solutions 

o 

of  the  following  matrix  integral  equations 


'“O 


N  T 

~Y  ^(t.z)  -f  ^  Kci(t,  u)  W,(u,z)du 


(u,z)du  =  K  (t,z) 
L 


N  Tf 

"T"  W«(fc»  z)  +  ^  K  (t,  u)  W  (u,  z)du  =  Kc  (t,  z) 
2  o  J  co  £  o 


(3.15) 


Ti 


where 


K_  (t,z)  =  K  (t,  z)  -f  K  (t,  z) 


and 


KcJ1*2)  =  K  (t,  z) 
o  c 


(3.16) 


These  equations  are  Fredholm  Integral  Equations  of  the  second 

kind,  and  can  be  solved  by  the  methods  given  in  34  and  35.  Following 

this  procedure,  the  optimum  decision  equation  in  terms  of  the  matrix 

filters  W,(t,  u)  and  W  (t,  u)  is 
I  o 


(t,  u)  r(u)  dt  du  - 


WQ(t,  u)  r(u)dtdu^  if 
(3.17) 


This  equation  can  be  rewritten  in  terms  of  a  single  matrix 


filter  M^(t,  u)  as 


J. 

L  =  ^  r  (t)  MA(t,  u)  r(u)  dtdu  <  Nq  n' 


(3.18) 


T. 

i 

where  the  integral  equation  for  MA(t,  u)  is 
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M^(z,  w)  Kn(w,  u)  dwdz  =  Kg(t,  a) 


(3.19) 


Equation  (3.10)  is  the  optimum  ''coupled  space-time"  decision 
equation  for  detection  of  random  signals  in  noise.  Interpretation  of 
equation  (3. 10)  yields  the  optimum  system  configuration,  the  following 
section  will  show  several  interpretations.  The  decision  equation  is 
written  in  terms  of  the  correlation  properties  of  the  incident  processes 
and  the  filtering  (space-time)  properties  of  the  array;  the  system 
requirement  will  reflect  this  same  dependence. 


3.  2.  System  Configuration 

The  previous  section  developed  the  optimum  decision  equation 


(3.18)  based  upon  the  maximum  likelihood  principle  for  binary  hypotheses. 


In  the  development  the  matrix  M^(t,  u)  was  referred  to  as  a  filter. 
Interpretation  of  equation  (3. 18)  yields  the  optimum  system  configurations. 


with  one  such  system  employing  the  M^(t,  u)  filter  to  be  discussed  first. 

Equation  (3.18)  can  be  rewritten  as 

T  T 

f  f 

L  =  ^  r  (t)  ^  M^(t,  u)  r(u)  dudt  (3.20) 

T.  T. 

i  i 


In  this  form  the  received  data  vector  r(t)  is  passed  through  a  lattice 
filter  bank  (convolution)  and  is  then  correlated  with  the  unfiltered  data 
producing  the  test  statistic  L.  The  matrix  filter  is  a  K-lattice  bank 


(for  K  receiving  elements)  of  filters.  The  system  configuration  is  shown 

in  Figure  3-1.  If  the  incident  processes  are  non- stationary,  the  filters 

37 

required  are  time  varying  and  can  be  realised  in  the  Kalman  form; 

* 

stationary  processes  at  the  input  would  require  linear  time  invariant 


filters. 


Another  system  configuration  for  the  optimum  receiver  can  be 


found  by  dividing  M^(t,  u)  into  two  filters  defined  by 


A  f 

x(t,  u)  =  \  Hw(t,z)  Qn(s,  u)dz 


(3.  21) 


where 


C  K  (t,  z)  Q  (z,  u)dz  =  6(t-u)I 
J  »  £  ~ 


(3.22) 


Substitution  of  (2.  31)  into  the  integral  equation  and  using 
equation  (3.  22)  fcr  the  inverse  noise  kernel  matrix  Qn(z,  u)  yields  the 
integral  equation  for  H^t,  z). 


i 

Kz(t,  u)  =  ^Kz(t,x)  +  Kn(t,x)J  Hw(x,  u)du  (3.23) 


34  35  14 

This  equation  is  the  KkK  matrix  Wiener-Hopf  equation,  ’  *  and  the 

matrix  filter  H^(t,  u)  is  the  minimum-means- square-estimation  filter 
for  the  K-vector  processes  z(t) 


z(t)  =  C  H  (t,  u)r(u)du 

—  mmse  j  w 


(3.24) 


structure,  and  of  the  steering. 


1-2.  Two-Combiner  Receiver  Structure,  interpretation  2. 


The  random  processes  are  said  to  be  stationary  if  the  observa 

tion  time  T  =  T^~  T.  is  greater  than  the  reciprocal  bandwidth  of  the 
26 

random  processes.  If  such  is  the  case,  Fourier  transformation  of 

equation  (3. 19)  yields  the  matrix  frequency  filter  Sq(u)  in  terms  of  the 

process  power  spectra  as 

S.  <uj)  =  S  (u,)S0(u>)  S  («,} 
z  4  n 

and 

Sq(iu)  =  Sz(m)  S^tu)  (3.27) 

where  S^(uj),  S^{ io),  S^(uj)  were  defined  previously. 

~  ~  14  34  35 

Employing  the  integrated  transform  *  ’  in  equation  (3. 18), 

and  using  Parseval's  theorem  gives  the  test  statistic  L  in  frequency 

form. 

CO 

L  =  C  RTtiu)  S  (ui)  R(a>)  ~  (3.  28) 

J  9.  -  ^TT 

-CO 

Substitution  of  equation  (3.  27)  into  (3.  28)  gives 

L  =  C  RT(u,)  s:\rn)  S  (w)  S  -1(uj)  R(oj)  ~  (3.29) 

j  J.  z_  n  -  ztt 

The  general  stationary  receiver  structure  is  given  by  equation  (3.  29).  . 
K  element  case.  A  simplified  structure  for  2  elements  is  shown  in 
figure  3-3, 

Paralleling  the  development  of  equation  'A-21)  through  (3.26) 


for  the  stationary  case,  gives 


where 


ZM-  S(»)S-lT<“lR<“»  (3-30) 

-  Z  1 

X(u>)  =  S“1(to)  R(j»)  (3.31) 


This  system  shown  in  figure  3-4  is  the  frequency  form  of  the 

structure  shown  in  figure  3-2.  It  will  be  noted  that  this  implementation 

again  requires  two  matrix  filters  and  combiners.  Simplification 

results  when  the  effects  of  array  steering  are  exploited  and  utilized 

14 

for  the  single  signal  case. 

It  was  shown  previously  that  array  steering  caused  the 
received  waveform  vector  to  undergo  phase  shifting  due  to  array 
steering  delays.  Thus,  after  steering,  the  signal  components  are 
identical  and 

S  .(id)  =  S  (<d)1  1T  (3.  32) 

z 1  z  —  — 

and 

V"»  ■  On'yH  ,3-  33> 


Thus,  equation  (3.  30)  becomes 


Fig.  3-4.  Stationary  Two-Combiner  Receiver  Structure,  interpretation  2. 
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32 


S2(w)_1T  Sj^uu)  =  T(uj)_1T  sjfa) 

S  (uj)1TI  =  T{uu)1T  S'ltiu)  S .» 
z  —  —  n  i 

but 

Sli(uJ>  =  sz»  +  Sn» 


Thus 


S  (io)lTI  =  T(u3)1TI+  T(uj)1T  S'huj)  S  .(id) 
z  n‘  z  ‘ 


After  simplification, 

Sz(u,)  =  [Sz(tu)  I  f  Sn«(u,)  *  X]  T((U) 

i  j 


Solving  for  the  scalar  frequency  filter  T(tu)  gives 


T(u>) 


sz(«,) 


K  K 
v’  r1 


i+s>>IX  s>> 

i  j 


S  (id) 
z 


1+  N(u>)  S  (uj) 
z 


Thus  the  most  easily  implemented  post- steering  receiver  structure 
passes  the  steered  data  through  the  KxK  matrix  bank  of  frequency 
filters  sums  the  outputs  of  the  filter  bank  tc>  form  the  scalar 

Y(id),  then  passes  Y(w)  through  the  scalar  linear  system  T (to) -  The 

a 

resulting  estimate  Z(i j?)  is  correlated  with  Y(oj)  to  form  the  test  sta- 
14 

tistic  L.  This  system  is  shown  in  figure  3-6  and  illustrates  the 
dependence  of  the  receiver  structure  on  incident  process  spatio- 
temporal  correlation  properties,  array  geometry,  and  array  steering. 


Mkrtx  Fitter* 


Figure  3-6.  Simplified  Single  Combiner  Receiver  Structure 


The  optimum  receiver  structure  can  be  found  in  the  above  manner  for 
a  given  incident  process  structure. 

The  above  sections  have  developed  the  optimum  decision 
equations  for  the  space-time  detection  of  random  signals  in  noise  and 
have  shown  the  optimum  system  configurations  that  result  from  an 
interpretation  of  these  decision  rules.  The  resulting  systems  were 
developed  for  the  non- stationary  and  stationary  steered  case  and  were 
shown  dependent  upon  the  space  time  properties  of  the  incident  signal 
and  noise  fields,  and  array  filtering  characteristics.  The  next  chapter 
investigates  the  performance  of  these  systems. 


CHAPTER  IV 


SYSTEM  PERFORMANCE 


Once  the  optimum  receiver  structures  have  been  found,  the 
performance  of  these  systems  must  be  evaluated.  A  performance 
measure  for  the  preceding  systems  will  be  established  and  the  nature 
of  the  test  statistic  L  will  be  investigated. 

To  determine  system  performance  characteristics,  one  needs 
the  probability  density  of  the  test  statistic  L,  under  each  hypothesis. 
From  this  density,  false  alarm  probability,  P_,  probability  of  detection, 

X 

Pp,  and  error  probabilities,  P^,,  can  be  defined  and  computed,  thus 

...  ,  ,  11-14,26 

specifying  system  performance. 

The  distribution  of  L  under  each  hypothesis  can  be  obtained  in 

two  ways.  The  first  is  a  direct,  exact  probability  density  calculation; 

the  second  uses  a  guassian  approximation  to  the  true  distribution  of  L. 

In  the  approximation  approach  the  density  of  L  is  taken  to  be  gaussian, 

with  true  mean  and  true  variance  computed  directly  from  the  decision 

equation.  The  validity  of  this  approximation  has  been  considered  and 

8,  36  36 

established  by  Price  and  Green 

In  both  cases,  the  density  of  the  test  statistic  under  each 
hypothesis,  and  resulting  performance  probabilities  will  be  shown  to 
be  dependent  upon  the  spatio-temporal  correlation  properties  of  the 
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i  b 

incident  random  processes,  the  array  geometry  and  steering,  and  the 
receiver  structure.  Analysis  of  the  performance  measures  will  yield 
guidelines  to  enable  one  to  improve  the  over-all  system  performance. 

In  this  chapter,  the  true  density  will  first  be  established,  then 
the  true  mean  and  variance  under  each  hypothesis  will  be  computed. 

Using  the  true  first  and  second  moments  the  gaussian  approximation 
will  be  made,  invoking  the  Central  Limit  Theorem;  and  from  the 
approximate  densities  the  system  performance  criteria  will  be  developed. 


4-1.  True  Distribution  of  the  Test  Statistic. 

The  optimum  space -time  decision  function  for  the  detection 
of  random  signals  was  given  by  equation  (3. 18)  as 


L  =  ^  r^(t)  M^(t,  u)  r(u)  dt  du 


(4. 1) 


T. 


Expanding  the  continuous  quadratic  form,  the  decision  equation  becomes 

Tf  K  K 

L  =  S  r^t)  MA„(t,  u)  rj(u)dtdu 

T.  1  J 

i 

K  Tf 


X  S  ri^  ri*u*  dt  du 


1  T. 
i 


K  K 

+  X  X  u)  rj(u)  dt  du 

i  j 


(4.2) 


Defining 


Z  (t)  =  V  MA..(t,  u)  r.(u)  du 
^  J  ij  J 
T. 


K  Tf  K  K  Tf 

L  =  y  ^  ri(t)  z^{t)  dt  +  y  y  ^  r.(t)z.(t)dt 

i-1  T.  i  j  T. 

i  ...  1 


Let  the  ij  integral  product  of  (4.  4)  be 


z^(t)  dt 


(4.3) 


(4.4) 


(4.5) 


Thus 


L 


K 


K  K 


i  i  j 

i*3 


L.. 

U 


=  Lll+L22+-'-  Lkk+L12  +  '”  Lk'lk 


(4.6) 


Examining  a  typical  term  in  (4.  6)  shows  that  each  L  term  is  an  infinite 
dimensional  inner  product.  This  is  shown  by  taking  the  general  term, 
Lmpj  and  expanding  it  in  a  series  of  orthonormal  basis  functions  using 
the  Karhunen  Loeve  expansion. 


L 

mp 


\  r  (t)  z  (t)  dt 
J  m  P 
T. 


(4.7) 


where  1  £  m,  p  £  k  and  the  processes  can  be  written 


2b 


In  the  limit  the  products  are  infinite  dimensional. 

In  order  to  find  the  density  of  L,  one  must  find  the.  density  of 
2 

the  sum  of  K  inner  products.  The  mathematic  is  quite  formidable, 

40  41 

as  the  products  Ln. .  are  non-gaussian  and  not  independent.  *  Thus 

27  38  39 

conventional  techniques  ’  *  are  not  applicable.  The  complete  density 

of  (4. 12)  is  not  investigated  here  but  in  order  to  illustrate  the  nature  of 
the  test  statistic,  the  present  section  will  obtain  the  exact  density  of  a 
gaussian  inner  product.  Such  a  density  is  not  without  utility,  since  it 
represents  the  test  statistic  density  computed  for  a  single  receiving 
element  system.  The  density  will  be  evaluated  in  terms  of  the  eigen¬ 
values  and  eigenfunctions  of  the  incident  processes  and  the  filtering  of 
the  frequency  filtering  sections.  The  Karhune  l-Loeve  expansion  will 
be  used  extensively. 


where  \r  is  the  ith  eigenvalue  of  the  process  r(t). 

:i 

The  inner  product  Ln  is 

Ln=rTz  (4.14) 

where  r,  £  are  gaussian  and  zero  mean,  as  in  the  process  model  taken 
in  Chapter  II.  Since  the  n-dimensional  vectors  are  expanded  using  the 


K-L  expansion,  the  co-ordinates  of  each  vector  are  independent,  while 


I 


there  is  correlation  between  the  elements  of  r  and  of  z.  This  is 


represented 


E[r.r.]  =  E[z.  z.]  =  0  Wj 

1  J  1  J  t  0  i  =  j 


E[rz.'J=0  it  j 

3  tO  i  =  j 


Using  this  property  of  independent  co-ordinates,  the  joint  density  function 
of  the  vectors  r,  z  is  given  by 


n 


where 


where 


and 


P(£»  z)  =  JT  p(r  »  z  ) 
i=l  1  1 


-1-  "l 


Ai^[cIxhCriZi3]  = 


E[r“]  E[r.z.] 


E[r.z.]  E[z^  J 


A?1  =  W. 

i  i 


(4. 15) 


_i  -V2  1  _ , 

p(r  ,  z  )  =  (2tt)  “  |  A.  |  exp(-j([r  z  ]  A-  [z.  ])  (4.16) 

11  1  1  1  L  J- 


(4.17) 


(4. 18) 


Thus 


The  joint  density  can  be  rewritten  as 


p(r,z)  =  ni  |  A.  |  2  exp -|(%.  rS  2W12.r.  z.  +  W22.  z.2) 


For  stationary  processes  A^  =  A.  all  i,  j 


n  n  n  n 

P (£■•*>  =  <^;>  Ul  2exP-|(WllIri2  +  2W12lro7'i+WE2lzi2> 

1=1  Si  i 


is  the  joint  density  of  the  r,  z  vectors.  The  exponential  summations 
can  be  redefined  as 

i  izi2= y2>  i  = Ln 


Thus 


p(r,z)  =  (2Tf)“n|A|  2exp(-|(WuX2  +  ZW^  £z  +  W ^  Y2)) 


Following  the  procedure  of  Miller,  generalized  spherical  co-ordinates 
for  the  r,  ^variables  are  introduced,  and  manipulated  to  give  the 
marginal  frequency  function  in  X,  Y,  cp  for  the  definition 


z_  =  X  Y  cos  cp  Oscp^TT 


giving 


F(X,  Y,  0)  = 


,VVWVY  •  ^  -^W11X  +W22Y)  "^XYcoscp 

(XY)  (XY  sincp)  e  ^  e 


zn“2|A|2  n|n)r(^)n|) 


(4.  20) 


Integrating  (4.  20)  over  the  range  of  the  new  variables 


JL.^  =  X  Y  cos  cp 


§  =  XY  sincp 


a  =  X 


0 


0  5a<® 


,  the  Jacobian,  gives 


h}  +  I2 

n  3 


(4.  21) 
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p(Ln)  =  ^  f(X,Y,cp)Jdad§ 


o  o 


With  K  (• )  the  modified  Bessel  function  of  the  second  Kind  and 

v 

order  v,  the  density  of  for  n  co-ordinates  in  the  expansion  is 


p(L  )  = 


n-1  -W.,L 

\h  |T“  e  12  n 

1  n' 


n 


n-1 

,1  .,2  ,  .  ,'2 


ITT  k-l<lLnlA^2>  (4-22> 


T(jn)Z  |  A  |  (Xu\22) 


The  matrix  A  is  found  from  (4. 17)  and  (4. 18)  in  terms  of  the 
space-time  correlation  properties  of  the  incident  processes,  array 
steering  and  filtering. 


ECriri 3  =  ^  Kr(t,  u)  ^(t)^(u)  dtdu 


T. 

l 


EO^]  =  E[z.r^]  =  ^  Kr(t,  v)  MA(v,  u)ijf.{t)  ^.(u)dudtdu 


T. 

i 


E[z^  ]  =  \^y\  Kr(v,  w)  MA(t,  v)  MA(w,  u)  ^(t)  \Jf.(u)  dvdwdtdu  (4.23) 
Ti 

In  terms  of  frequency  spectra  using  integrated  transforms 

•  ■»  , 
dw 


E[r.  ]  =  Sr(ui)  |\lf.(ju))|  ~ 


E[r.z.]  =  E[z.r.] 


=  ^  Sr(u))MA(iu)|i!f.(jw)| 

~  03 
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«  2  2 

^  Sr(tu)  }MA(io)j  I$.(j(0)l 

-00 

where  ^ 

^.(ju))  =  ^  \!f.(t)  e  ^  dt 

T. 

i 

Thus  the  density  is  dependent  upon  the  space-time  structure  of 
the  incident  signal  and  noise  processes  in  the  above  manner.  The  density 
given  in  (4.  21)  with  (4.  23)  can  be  written  under  each  hypothesis  by  appro¬ 
priate  substitution  of  (2. 19),  {2.  20)  into  the  correlation  function  in 
equations  (4.23). 

In  summary,  the  continuous  density  of  L  under  either  hypothesis 

is  found  by  obtaining  the  limit  of  the  density  of  as  n  grows  without 

bound.  An  approximation  can  be  realized  by  trancating  the  infinite  series 

at  n  and  using  n  co-ordinates.  The  density  given  above  is  the  density  of 

2 

a  single  inner  product.  The  true  density  must  bp  found  for  the  sum  of  K 
such  products.  Since  each  of  the  inner  products  is  non-gaussian  and  not 
independent,  the  analytical  effort  is  not  very  tractable  and  has  not  been 
investigated  here.  The  density  of  the  single  inner  product  shows  the 
nature  of  the  test  statistic  L.  and  in  a  following  section  the  sum  of  many 
such  products  is  allowed  to  converge  to  the  gaussian  distribution,  thus 
allowing  system  performance  to  be  expressed  analytically.  The  next 
section  of  this  paper  obtains  the  true  mean  and  variance  of  the  complete 


continuous  test  statistic. 
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4-2.  Calculation  of  True  Mean,  and  True  Variance  of  L 

The  true  first  and  true  second  order  moments  of  L  (mean  and 
variance)  under  each  hypothesis  can  be  found  by  direct  application  of 

t  ■ 

the  expectation  operator  to  the  decision  equation.  Under  the  received 
vector  r(t)  is  taken  to  be 

r(t)  =  z(t)  +  n(t) 

Thus  taking  the  expected  value,  the  mean  of  equation  (3.10)  under  is 

Tf  _ 

E[L|Hj]  =  e[  CC  rf(t)  MA(t,u)  r(u)  dtdu]  (4.24) 

T1 

and  since  the  integrand  of  (4.  24)  is  a  Hermitian  quadratic  form  of  the 

T  40 

form  R  QR,  the  form  can  be  rewritten  as 


RTQR  =  Tr[QRRT] 


(4.25) 


where  T^  iB  the  trace  operator. 

Thus  (4.24)  becomes,  after  commuting  the  expectation. 


i.  _ 

ECLjHj]  =  Tr  g  MA(t,u)  Efc't)r(u)]  dtdu 


i 

=  Tr  ^  M^(t,  u)  K^(t,u)dtdu 


(4.26) 


K  (t,u)  =  K  (t,u)  + K  (t,u) 

x  Zi  n 


with 
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In  a  similar  manner  the  expected  value  of  L  under  Hq  is  found 


A 

E[LjHQ]  =  Tr  MA(t,  u)  Kn(t,  u)  dt  du 


(4.  27) 


The  variance  of  the  text  statistic  is  found  from  the  definition 


of  the  variance. 


VarCLlH^  =  E[L21H.]  -  E2[L|H.] 


(4.28) 


Squaring  (3.18)  to  find  the  mean- square  of  L  under  H. 


1»2  =  r^(t)  MA(t,  u)  r(u)  r^(x)  M^xf  y)  dt  du  dx  dy  (4.  29) 

i 


and  applying  property  (4. .25)  and  the  property  that  matrix  products 

j  40 

commute  under  the  trace  operator 

Tr  (ABC)  =  Tr  [any  order] 


gives 


L  = 


=  •••  Trj^MA(t,  u)M^x,  y)[r(t)r  (u)r(x)r  (y)]  J  dxdy  dt  du 


Applying  the  expectation  under  IT  and  the  expectation  theorem  for  the 

33  26 

product  of  gaussian  variables  }  it  follows  that 


E[L  |H  ]  =  2Tr 


MA(t,  u)  K^u,  x)  Ma(x,  y)  K{(y,  t)  dt  du  dx  dy 


(4.  30) 


+  Tr 


MA(t,  uJK^t,  u)Ma(x,  y)Ki(x,  y)  dt  dudxdy 


The  second  set  of  integrals  in  H..  30)  is  the  E  [L|Hj]  which, 
when  taken  with  (4.  29)  gives 


T. 

i 


MA(t,  u)  K.(u,  x)Ma(x,  y)K;(y,  t)dtdudxdy 

(4.  31) 


Var[Lj!T]  =  2Tr  "*&*'’*  “/  *•*£'**»  , — £  -  l 


Thus  the  true  mean  and  true  variance  under  each  hypothesis  are 
given  by  equations  (4.27,  31).  These  equations  give,  in  general,  the 
first  and  second  moments  for  both  stationary  and  non- stationary  processes 
If  (r(t)}  is  a  stationary  process,  the  equations  (4.  27,  31)  can  be  written, 


14  34  35  14 

using  integrated  transforms,  *  ’  and  Parseval's  theorem*  in  terms 


of  the  process  power  spectra. 


E[L|H.]  =  [  S^oj)  Sz(u>)  S’\«)  S.(u>)  (4.32) 

and 

03 

Var  [L|H.]  =  2Tr  ^  |  S'V)  S^w)  S^(u,)  S.(»j»)  j 2  ~  (4.33) 

-00 

The  above  equations  for  the  true  mean  and  variance  of  L  under 
each  hypothesis  will  be  used  extensively  in  the  gaussian  apprlximation 
to  the  true  density  that  follows  and  in  evaluating  system  performance. 


4-3.  Gaussian  Approximation  for  the  Test  Statistic  Distribution 

The  density  of  the  single  element  test  statistic,  L  under  either 
hypothesis  was  given  by  equation  (4.  22)  As  can  be  seen,  this  density  is 
quite  complicated  and  difficult  to  manipulate  analytically  to  obtain  exact 


system  performance.  Such  difficulty  can  be  overcome,  however,  by  the 


4? 


the  use  of  an  approximation  to  the  exact  density.  One  such  approximation 
is  a  truncation  of  the  infinite  series  in  (4.  22)  and  the  use  of  numerical 
techniques  to  obtain  system  performance.  Another  approach  is  a  direct 
approximation  of  the  true  density  by  substitution  of  a  gaussian  density 
for  (4.22).  In  this  approximation  L  is  taken  to  be  a  gaussian  variable 
with  mean  given  by  the  true  mean  of  L  (computed  in  (4.  27,  28)  and  with 
variance  given  by  the  true  variance  of  L  (computed  in  (4.  31)). 

Justification  for  this  approximation  can  be  seen  by  noting  that 

the  test  statistic  L  was  found  in  IV-2  to  be  a  kth  order  inner  product, 

,  ,  .  ,  .  ,  ..2  .  2  .  ,  .  27, 40, 41 

and  as  such  is  the  summation  of  K  non-gaussian  X  variables. 

It  is  known  that  the  density  of  a  summation  of  non-gaussian  inch  variables 

converges  to  the  gaussian  distribution  quite  rapidly  through  the  proper¬ 


ties  of  the  Central  Limit  Theorem. 


27,38 


This  convergence  is  most 
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rapid  in  the  vicinity  of  the  mean  of  the  density.  Price  ’  ,  Green  , 

42 

and  Middleton  have  investigated  the  gaussian  approximation  quite 
extensively  and  ha‘/e  shown  its  validity  in  computing  system  performance. 

The  gaussian  approximation  for  the  test  statistic  has  been  inves¬ 
tigated  for  the  "threshold"  or  coherently  undetectable  case  by 
42  36 

Middleton  and  Price  .  Their  studies  have  shown  that  the  approxima¬ 


tion  yields  useful  results  for  cases  in  which  the  signal-to- noise  ratio 
(SNR)  is  low.  If  a  large  average,  signal  power  is  available,  implying 
high  SNR,  system  detection  requirements  are  usually  not  critical. 


However,  for  low  SNR  the  effects  of  array  and  incident  process  spatio- 
temporal  correlation  properties  become  important,  since  it  is  usually 
necessary  to  attain  nearly  optimum  system  configuration  and  performance, 
in  order  to  meat  system  detection  capability  specifications.  The 
gaussian  approximation  furnishes  performance  analysis  methods  for  such 
cases  of  low  SNR.  Thus  the  approximation  is  most  valid  where  it  is 
of  greatest  utility.  Helstrom  has  examined  another  approximation  to 
(4.  22)  using  the  Gram-Charlier  series  expansion  and  found  comparable 
results. 

Thus  the  density  in  (4.  22)  can  be  approximated  as  gaussian 
with  mean  given  by  {4.  27,  28)  and  variance  by  (4.  31).  The  following 
system  performance  criteria  analysis  and  subsequent  case  studies 
employ  the  gaussian  approximation  of  the  test  statistics. 

With  the  density  of  L  gaussian,  the  error  probabilities  of  the 
system  can  be  defined  by 

00 

prob.  false  alarm  =  =  V  p(L|HQ)dL 

Y 

=  “Kt1  +  t) 

0 


where  y  is  the  threshold,  and  similarly 


prob.  of  miss  =  P. 


=  f  P(L|Hj) 


=  erf  { 


■»  ° 
Inn  , 


{4.  35) 


where  erf(x)  =— — 


$  exP<:T: 


erfc{x)  =  1  -  erf(x) 


P  ,  the  probability  of  detection,  is  defined  as 

p  =  |~  p  =  erfc(^LU--^) 

D  1  M  '  6  2  1 

o 

The  parameter  6q  is  related  to  the  output  signal  to  noise  ratio 


So  =  SNKo 


where 


(E[LjH  ]  -  E[L|H  ])  2 

(SNR  f  = - - - —  1,  =  6 

°  (VarCLlH^VartLl^])'2 


(4.  36) 


The  validity  of  this  measure  of  system  performance  has  been 
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found  by  Price  .  The  numerator  of  the  6q  parameter  is  a  measure  of 
the  distance  between  the  means  of  the  densities  under  each  hypothesis, 
while  the  denominator  is  a  measure  of  the  variance  or  spread  of  the 
densities.  Thus,  system  performance  will  be  enhanced  by  increasing 


the  distance  between  the  means,  and  degraded  by  an  increasing  variance. 


The  6q  parameter  will  be  irsed  as  measure  of  system  performance 
in  the  remainder  of  this  paper. 

Using  equation  (4.  36)  for  6  2  and  equations  (4.  27),  (4.  28), 

v  O 

(4.  31)  for  the  parameters  in  (4.  36),  the  performance  characteristics 

can;  be  written  in  general  as^ 

Tf  2 

tTr  S  MA(t'  Kz(t*  dr  du] 

62=  Ti 

o 


i 

'C^TlS5  MA(t*  *VU’ M^x*  y)Kx(y* l) dt  dudxdy] 

x  CTr  MA(t,  u)Kq(u,  x)Ma(x,  y)K^y,  t)  dt  dudxdyf) 


(4.  37) 
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The  parameter  6q  reflects  the  dependence  of  system  performance 
upon  spatio-temporal  structure  of  the  incident  processes,  array  geometry 
and  steering,  signal-to-noise  power  ratio,  and  temporal  filtering  and 
combining  techniques.  Equation  (4.  37)  is  the  general  expression  for  the 


receiver  performance  for  both  stationary  and  non- stationary  processes. 

If  the  incident  processes  are  stationary  as  defined  previously,  new 

2 

insight  into  the  nature  of  is  found. 

With  stationarity  and  binary  hypothesis  6.°  can  be  written  in 
terms  of  the  power  spectra  of  the  incident  processes,  filtering  require¬ 
ments  and  steering  delays  as 


and 


s«b! 


Si 


S  (uj)  =  S  1{U))  S  (oj)  S  ^od) 

U  1  z  n 


«  S1(u»)  s_(w)  Sr.1(w) 

n  z  n 


Using  these  relations  in  (4.38)  and  noting  that  Var £L|  H^]  SJ  Var[L.|H^J 
gives 

2  ^ECLjH^-ECLlH^}2 


6  =T 
o  2 


Var[LjHoJ 


'09 

4ti  \  s  (tu)  S  (id)  S~1(uj)  S  (uj)  |r 
2  j  n  z  n  z_  £TT 


for  a  steered  array  S  (uj)  =  S  (uj)  is  given  by  (3.  32)  and  S  (uj)  -  S  ,(uj) 

z  z*  n  ix 

by  (3.  33).  Thus 


and 


5  = 


^rSs>»iiTs;iiiTs>g 


s>)IIT  =  [l  s>>] 


K 


Tr[(Sii»2iT)2]  =  (H  S>0 


K  K 


i  i 


(4.  39) 


(4. 40) 


Substitution  of  (4.  40)  into  4.  39)  and  taking  the  trace  operator  inside 
the  integral  gives 


6o  •  1 (s>>  It  s>D  i 

i=ij=i 


(4.41) 
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Squat  ion  (4.41)  is  one  final  result  for  the  system  performance 

2 

index  5^ .  The  index  is  expressed  in  terms  of  the  input  signal  spectrum 

S  {(u),  which  may'nave  a  multipath  representation,  and  the  steered  noise 

matrix.  The  latter  contains  the  effects  of  array  geometry,  array  steering, 

and  spatio-temporal  correlation  properties  of  the  incident  processes  as 

well  as  system  interference.  The  elements  S_i  (id)  can  be  found  from 

ij 

2 

equation  (2.18).  Using  these  results,  one  can  calculate  the  index  6q  for 

various  values  of  array  steering  and  geometry,  as  well  as  different  space- 

2 

time  properties  of  both  the  signal  and  noise  fields.  Using  6q  one  can 
then  obtain  the  receiver-operating-characteristics  (ROC)  curves  to 
specify  system  performance  in  terms  of  error  probabilities.  Using  these 
curves,  the  designer  can  establish  the  effects  upon  system  performance 
of  various  parameters  in  the  system,  such  as  the  effects  of  the  number 
of  array  elements,  various  array  geometries,  multipath  structure,  and 
of  the  cross  element  space-time  correlation  properties  which  are  related 
to  the  array  steering  and  incident  process  structure  by  equations  (2.13) 
and  (2. 18). 

The  previous  chapter  developed  the  optimum  decision  equations 
for  the  array  detection  of  random  signals  in.  noise.  Interpretation 
yielded  a  number  of  system  configurations  (both  stationary  and  time- 
varying).  The  performance  of  these  optimum  systems  was  then  investi¬ 
gated  for  both  the  general  case  and  the  single  signal  (steered  array) 
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stationary  process  case.  Introduction  of  the  threshold  assumption 
allowed  further  simplification  of  the  performance  index  and  interpre¬ 
tation  to  illustrate  the  effects  of  steering  and  cross  link  properties  of 
the  random  fields.  The  following  section  will  investigate,  as  case 
st tidies,  the  performance  of  the  optimum  systems  to  determine  the 
effects  of  array  geometry  (element  spacing  and  position),  the  effects  of 
multipath  in  the  incident  signal  processes,  and  the  effects  of  cross¬ 
element  correlation  properties  of  the  random  fields.  Analysis  and 
interpretation  of  the  resulting  performance  curves  (ROC)  will  establish 
guidelines  as  to  techniques  for  enhancing  system  performance,  by  valid 
physical  consideration. 


CHAPTER  V 


CASE  STUDY 

Using  the  results  of  Chapters  HI,  TV,  and  the  model  of 
Chapter  II,  the  following  chapter  analyzes  in  detail  the  performance 
of  the  previously  derived  system  configurations  for  various  cases. 

The  first  case  study  will  investigate  the  detection  performance 
of  the  single  element  optimum  receiving  system  configuration.  The 
incident  processes  will  be  assumed  gaussian,  with  the  signal  portion 
of  the  received  process  consisting  of  a  two  path,  correiated-multipath 
random  process  with  "Markovian"  correlation  function.  The  system 
performance  is  evaluated  in  terms  of  the  effects  produced  by  the  delay 
between  the  multipaths,  and  the  division  of  signal  power  between  the 
two  resolvable  paths,  thus  establishing  the  effects  of  the  multipath 
structure  upon  single  receiving  element  system  performance. 

The  second  case  study  will  employ  a  two  element  receiving 
array.  In  this  study  the  effects  of  multipath  structure  upon  system 
performance  will  again  be  established.  In  addition,  the  effects  of 
element  geometry  and  cross-link  spatio  temporal  correlation  properties 
of  the  incident  random  processes  will  be  pinpointed.  The  physical 
aspects,  multipath  effects,  element  geometry,  and  noise  correlation 
properties  will  then  be  combined,  illustrating  a  method  for  enhancing 
system  performance  by  proper  selection  of  array  geometry. 
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Finally,  using  the  "best  1  geometry  found  in  the  second  case 
study,  the  effects  of  the  number  of  array  elements,  upon  performance, 
wilXbe  illustrated. 

Thus,  these  case  studies  establish  the  performance  character¬ 
istics  of  single  and  multiple  element  receiver  systems.  They  illustrate 
variations  in  system  performance  due  to  multipath  structure  in  the 
random  signals^array  geometry,  array  steering,  and  space-time 
correlation  properties  of  the  incident  random  processes.  In  general, 
the  case  studies  establish  the  role  of  spatio-temporal  correlation 
properties  of  the  incident  processes,  in  terms  of  coupling  by  the  array 
receiving  system,  for  the  random,  multipath  signal  and  noise  environ¬ 
ment. 

5-1.  Homogeneous  Interference  Assumption 

The  steered  noise  spectral  matrix  was  given  by  equation  (3.  33). 
If  the  noise  spectrum  at  each  receiving  element  is  assumed  to  be  the 
same,  the  noise  interference  field  is  said  to  be  a  homogeneous  field. 

Such  a  space-time  structure  of  the  noise  field  is  frequently  encountered 
in  both  the  sonar  and  radar  cases,  as  well  as  in  seismology. ’^The 
homogeneity  of  the  noise  is  represented  mathematically  by 

S  (uj)  =  S  (u>)  all  i,j 
n  ij  n  ij 

In  addition,  it  is  assumed  that  there  is  some  common  factorable 
noise  spectral  property  in  the  elements  of  the  noise  spectral  matrix. 


Thus,  equation  (3.  33)  can  be  written  as 


S,(u>)  =  S  (uj)  §  (tu) 
n*  n  - 


(5.1) 


where  the  §  (m)  matrix  has  elements 


G).j(uj) 


(5.2) 


Substitution  of  the  above  equations  into  the  threshold  detection  performance 


index,  (3.  39),  gives  the  new  index 


so 


■  sc$x£f  ."•<»)'  s 

n  i=lj=l 


-OO 


(5.3) 


Defining 


Sz(w)  A 

o  v  .  =  F  (uo)  =  signal-noise-ratio-function 
S  (tu)  o 
n 


and  the  array  gain  function 
K  K 


n  cp  ^  (uj)  =  G  (a))  =  array  gain  function 

i=lj=l 


one  can  write  the  performance  index  as 


=  C  r£(.)GZ(,| 

o  J  O  o  2tt 


(5.4) 


(5.  5) 


It  is  seen  from  equation  (5.  5)  that  the  performance  is  a  function  of  the 
input  signal-to-noise  ratio  and  an  array  gain  function  which  contains 


•oC 


the  enhancement  due  to  the  steering,  frequency  filtering,  and  combining 


as  well  as  the  coupling  produced  by  the  array.  It  is  the  index  in  (5.  5) 
that  is  employed  to  analyze  the  following  case  studies. 

5-2.  Case  I  -  Single  Element  in  Random  Multipath  Environment  -  Single 


Plane  Wave  Coherently  Undetectable  Case 
For  the  first  case  study,  the  receiver  consists  of  a  single 
receiving  element.  The  system  decision  equation  (3.18)  for  detection 


of  the  random  stationary  signal  z(t)  becomes 


co 

L  =  ^  R*(uj)  Sq(uj)  R{u>) 


(5. 6) 


with  the  filter  Sq(uj)  given  by 


sQto)  = 


Sz{m) 


S  (oj)CS  \uj)  +  S  (u>) 
n  z  n 


The  portion  of  Sq(ou),  [Sz(ui)/  S^(w)\  is  recognized  as  thenWiener'test5.mation 


,...  43, 14 

filter. 


This  system  configuration  is  shown  in  the  figure,  5-1,  with  the 


signal  spectrum  S  (u>)  being  found  below.  The  performance  of  the  system 
z 

is  found  from  (5.  3),  where  the  array  gain  equals  1  and  S  ,(uj)  =  S  (iu). 

n'  n 


The  performance  index  is 


S  {(!)). 


6°  =  S  C  S  (u>)) 


(5.7) 


.  .V.v.kww*. 
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where  uj  is  the  3dB  bandwidth  of  the  random  signal  process  and  P  is 
c  z 

the  average  power-  Thus  S  (id)  in  (5.  6)  and(5.  7)  is 

z 


2  2  2uj  P 

S  (id)  =  (A  +  B  -j-2ABcosujA)  — ~ — 

Z  L  c. 

OJ  T  U) 

c 


(5.11) 


The  system  performance  index  is 


C  /"A  +B  +  2AB  cos  ojA)  2ujcR-\  duu 

=  - r— - —  )  a?  (5-12) 


2  2 
U)  +  U3 
C 


A  plot  of  the  effective  signal-noise-ratio,  6q/SNR,  versus  the 
normalized  %  of  power  in  the  second  path  is  shown  in  figure  5-2  for 
various  values  of  multipath  delay. 

The  plot  shows  that  performance 

1.  degrades  as  the  multipath  delay.  A,  increases 

2.  degrades  as  a  larger  portion  of  the  signal  power  is 
transmitted  via  the  second  path. 

The  limiting  curve  of  "worst  performance"  is  the  curve  A-*®, 

while  the  "best  performance"  is  for  A=0,  implying  no  multipath  structure. 

The  physical  effect  of  the  multipath  structure  is  to  distribute 

the  signal  energy  over  a  larger  amount  of  time.  This  can  be  seen  by 

assuming  that  during  the  signal  band,  Tf-  T^  =  T.  for  no  multipath,  E, 

a  fixed  amount  of  energy  is  transmitted.  The  signal  energy  per  time 

E 

unit,  available  to  the  receiver  is  —  and  the  input  signal-to.-noise  ratio, 
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Figure  5-2 


with  N^  representing  the  noise  power  is 

SNRi=  ilr 

P 

If  two  paths  exist  with  A%  of  the  energy  transmitted  via  the  first  path 

and  B%  transmitted  via  path  two  at  time  A,  the  total  energy  available  is 

still  E,  since  A-j-  B  =  1,  but  the  time  interval  of  total  transmission  is  now 

T-fA.  The  input  signal-to-noise  ratio  is  lower,  shown  by 

AE-fBE  _  E 
(T+A)N  TN 

The  effective  SNR  given  by  the  performance  index  (5. 12)  shows 
quantitatively  the  above  effects.  In  figure  5-2,  the  index  decreases  due 
to  increasing  A,  shown  qualitatively  above,  and  is  "worst"  for  A  -*■<». 
Best  performance  is  achieved  when  A  is  small;  the  signal  energy  con¬ 
centration  being  greatest  in  the  shorter  transmission  or  observation 
time  period. 


The  degradation  due  to  more  equal  distribution  of  energy  over 
the  two  paths  with  constant  A  is  seen  by  examination  of  the  SNRj  for  a 
single  path.  For  A  %  of  E  transmitted  i.n  the  first  path  during  T,  and  B% 
in  the  second,  the  SNR  for  each  path  is 
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The  receiver  in  effect  must  process  two  weaker  signals  in  the 
presence  of  constant  noise  power,  rather  than  a  single  strong  signal. 
Thus,  the  effective  SNR  for  each  path  is  lower  than  for  a  single  path, 
and  the  resulting  performance  degrades.  This  effect  is  shown  in  the 
figure  by  the  decreasing  performance  index  for  increasing  power  ratio 
(or  increasing  energy  in  the  second  path).  Thus,  the  performance  for  a 
single  element  receiver  is  optimum  when  the  multipath  structure  is 
such  chat  a  majority  of  the  available  signal  energy  is  transmitted  in  a 
single  path.  If  the  power  does  split  appreciably,  the  best  performance 
occurs  for  small  delay  A. 

The  performance  for  the  single  element  receiver  and  given 
environment  cannot  be  improved  without  some  influence  upon  the  multi- 
path  structure  of  the  information  channel.  As  this  is  almost  always. not 
the  case,  since  the  multipath  structure  is  normally  a  natural  phenomenon 
of  the  channel  and  not  under  the  designer’s  influence,  the  performance 
of  a  single  element  is  almost  entirely  fixed  by  the  channel,  with  little 
improvement  capability  available.  Figure  5-3  shows  ROC  curves  for 
the  single  receiver  element  normalized  SNR,  and  various  multipath 
structures,  illustrating  the  performance  variations  noted  above. 

5-3.  Case  II  -  Two  Element  Array,  No  Multipath  -  Coherently 
Undetectable  Case 

Turning  now  to  the  array  detection  case,  the  following  case 


study  analyzes  the  performance  of  a  two  element  array  where  the 
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random  signal  has  no  multipath  structure.  In  this  situation  the  steered 
signal  spectrum  matrix  (3.  32)  is 

S  .(«,)  *  S  (uj)  11T 

Z  *  Z 

The  steered  noise  spectrum  matrix  is  (3.  33) 

S>>  -  tSnO 

The  elements  Sai..{tu)  are  given  by  equation  (2.  8)  and  the  waveform  delay 
is  given  by  (2.13).  Thus 

r 


s>>  ■ 


s„  (ul) 


n 


H 


S  (uj)  e 
“12 


-jtUT 


S  (uj)  eJU)T  S  (uj) 
n21  n22 


(5.13) 


Application  of  the  homogeneous  noise  field  assumption,  and 
assuming  that  the  correlation  properties  of  the  unsteered  noise  matrix 
are  such  that 


Sn-» 

y 


=  S  » 


N 


i  =  3 


(5.14) 


=  y  +  sc(w)  i  i  3 


where  Sc(m)  is  the  common  link  power  spectrum,  the  steered  noise 
matrix  is 


<5.15) 


6? 


where  A  is  defined 


A  Sc<“> 
A  = 


N 

—  +  S  (m) 

L  C 


(5.16) 


The  optimum  system  configuration  in  single  combiner  form  is  shown  in 
figure  5-4. 

Putting  equation  (5. 15)  into  the  form  of  (5. 1)  gives  the  common 


spectral  occupancy  S^cu)  and  the  §(iu)  matrix. 

Sn{w)  =  T  +  Sc{U,) 


(5.17) 


cp*(u))  = 


Finding  cp*^(m) 


-JU)T 


1  Ae 


A  eJU,T  1 


co’^uu)  = 


-Ae' 


-Ae“JU,T  1 


JU>T 


(5.18) 


(5.19) 


1  -  A 


and  the  array  gain  function  (eq.  5s- 4)  becomes 


G„<«)  = 


2(1  -  A  cos  uj  t) 
l-A2 


(5.20) 


For  mathematical  simplicity,  assume  that  the  noise  is  entirely 
N 

white,  with  S  (to)  =  ;  thus 

O  d* 


A  =  i-  ,  S  (u>)  =  N 
2  n  o 
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Substitution  of  these  terms  into  the  performance  index  equation 
(5.  3)  and  taking 


S2(ui)  = 


2ui  P 

C  2 

2  2 
0)  +  U) 

c 


in  Case  Study  I  gives 


2  »  .2«-cP-'(<«c2+»2)2  l 

o  "  )  v  N  Jl!“2 


COS  U)T 


v 

y  2tt 


(5.  21) 


Integration  of  (5.21)  and  plotting  the  error  probabilities  for 
various  values  of  t  gives  the  ROC  operating  curves  given  in  figure  5-5. 
These  curves  show  the  variation  o'f  the  system  performance  with  changing 
element  separation,  related  by  equation  (2. 13).  The  "best"  performance 
is  realized  for  infinite  separation,  while  "worst"  performance  results 
for  t  =  0,  (no  separation). 

The  improvement  in  performance,  obtained  by  making  T  as 

large  as  possible,  can  be  interpreted  in  terms  of  the  "detectability 

resolution"  of  the  entire  space  time  receiving  system.  In  antenna  design 

(spatial  detection),  it  is  often  desirable  to  make  the  main  lobe  of  the 

beam  pattern  as  narrow  as  possible  while  simultaneously  reducing  the 

4  5 

side  lobes  (minor  lobes)  of  the  antenna  array  pattern.  *  These 
requirements  are  mutually  exclusive,  to  a  certain  degree,  since  narrow¬ 
ing  the  beam  usually  results  in  increased  sidelobe  level.  The  problem 
is  usually  solved  by  employing  an  array  of  elements  spaced  and  steered 
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Figure  5-5 


to  give  a  narrow  directional  main  beam.  The  power  received  by  each 

2  4  5 

element  is  then  weighted  (shaded)  to  give  a  desired  sidelobe  level.  *  * 

There  are  numerous  shading  schemes  to  obtain  certain  sidelobe  charac- 
2  44  4  5 

teristics.  3  *  3  In  these  methods  the  physics  (space-time)  structure) 

of  the  noise  field  is  not  takep  into  the  problem  analysis,  thus  the  resulting 

resolution  (detection  enhancement)  is  due  to  spatial  resolution  only.  It 

45 

has  been  found,  particularly  in  the  area  of  radio  astronomy  ,  that  '‘best" 
resolution  is  achieved  by  widely  spacing  the  elements  of  the  array.  As 
the  element  spacing  increases,  ihe  resolution  increases  markedly. 
Multi-element  arrays  with  large  element  separation,  such  as  in  inter¬ 
ferometers,  are  widely  employed  in  much  of  present  day  radio  astronomy 
research.45 

In  the  case  of  optimum  space-time  detection,  the  analysis 
includes  the  effects  upon  detection  capability  that  arise  from  considera¬ 
tion  of  the  space -time  noise  fields  as  well  as  the  spatial  array  beam 
pattern.  Here  also  some  of  the  same  characteristics  of  the  receiving 
systems  are  wanted:  narrow  directional  main  lobe,  low  sidelobe  level, 
and  effective  noise  field  rejection.  The  latter  requirement  is  added 
because  the  array,  together  with  the  steering,  combiner,  and  lattice 
frequency  filtering  stages,  produces  effective  noise  rejection  and  signal 
enhancement  by  adjustment  of  the  sidelobe  levels.  The  noise  rejection 
and  signal  enhancement  result  from  the  optimal  system  configuration 
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des-eloped  in  Chapter  HI  and  discussed  there.  Thus  the  space-time 
analysis  yields  the  system  detection  performance  in  terms  of  both  the 
spatial  resolution  and  the  filtering  (frequency)  resolution  of  the  incident 
processes.  The  combined  effects  and  resulting  performance  index  can 
be  termed  the  total  array  "detectability  resolution,”  combining  the  disci¬ 
plines  of  antenna  resolution  theory  and  statistical  communications  and 
optimal  filtering  theory. 

The  increased  performance  capability  observed  in  figure  5-5 
with  increased  separation  of  the  elements  is  thus  the  result  of  sidelobe 
1)  noise  rejection  due  to  optimal  lattice  combining  in  the  frequency  filter¬ 
ing  stages;  2)  narrowing  the  main  lobe  (yielding  improved  3patial 
resolution)  due  to  increased  separation  of  the  elements.  The  element 
geometry  (separation)  and  steering  produce  a  spatial  beam  pattern  with 
a  narrow  directional  main  lobe;  the  optimum  frequency  filtering  and 
combining  stages,  dictated  by  the  space -time  structure  of  the  incident 
processes  and  the  array  geometry,  perform  sidelobe  rejection  and 
shading,  as  well  as  statistically  filtering  the  received  data.  Increasing 
the  separation  produces  a  narrower,  more  directional  spatial  beam 
pattern,  thereby  improving  resolution  in  space,  and,  since  the  new 
optimal  frequency  filtering  and  combining  stages  contain  the  new  spatial 
variables  as  parameters,  these  stages  produce  a  new  optimal  noise 
rejection  and  sidelobe  level  adjustment  pattern,  effectively  yielding 
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an  improved  detection  capability.  The  improvement  in  performance  due 
to  element  separation  can  be  seen  by  plotting  the  beam  pattern  of  the 
above  array  receiving  system.  In  figure  5-6  the  beam  pattern  for  a 
two  receiving  element  system  is  shown  for  the  above  interference  field 
and  two  values  of  element  spacing.  The  solid  curve  represents  a  spacing 
of'X/2,  while  the  broken  curve  is  the  beam  pattern  for  a  spacing  of  1.  OX.. 

As  the  spacing  is  increased  from  X./2  to  1.  OX.,  the  beam  pattern  changes 
in  two  ways;  the  main  lobe  narrows  and  the  side  lobe  level  increses. 

The  frequency  filtering  and  combining  stages  will  compensate  for  the 
increased  sidelobe  level  by  producing  optimal  noise  rejection- signal 
enhancement  in  frequency,  while  the  narrow  main  beam  will  produce 
improved  signal  resolution  in  space.  The  over -all  result  is  improved 
detection  performance  capability.  The  same  results  are  seen  in 
figure  5-7.  Here  the  beam  pattern  for  the  two  element  system  is 
plotted  for  two  spacings  and  a  signal  incident  at  30°.  The  increased  spacing 
gives  a  narrower  main  lobe  and  improved  detectability  resolution. 

Therefore,  consideration  of  the  detection  problem  as  a  combined, 
integrated,  space-time  problem,  rather  than  dichotomizing  the  problem 
into  a  spatial  (antenna)  problem  and  data  processing  problem,  yields  a 
more  concise,  broader  interpretation  of  design  requirements.  The 
above  case  study  shows  that  improved  system  performance  is  realized 
for  a  two  element  array  in  a  non-multipath  environment  by  making  the 
separation  of  the  receiving  elements  as  large  as  is  physically  possible. 


Coherently  Undetectable  Case 


The  preceding  case  study  has  shown  the  system  configuration 
and  performance  characteristics  for  a  two  element  array  operating  in  an 
environment  in  which  there  was  ho  multipath.  The  following  case  study 
will  investigate  the  system  performance  for  the  same  system  imbedded 
in  the  same  noise  field  environment,  but  with  the  random  signal  now 
having  the  multipath  structure  that  was  assumed  in  Case  I. 

The  system  configuration  is  the  same  as  that  in  figure  5-4,  but 

with  S  (10)  now  given  by  equation  (5. 11)  of  Case  I. 
z 

The  system  performance,  a  combination  of  input  process  SNR 
given  by  (5. 11)  and  array  gain  given  by  (5. 10),  is  given  by  equations 
(5.  3)  and  (5.  4)  as 
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(5.22) 


where  the  parameters  to  ,  P  ,  T  and  A  were  defined  previously. 

c  z 

Integration  of  equation  (5.  22)  and  numerical  evaluation  with  t 
fixed  and  various  values  of  A,  B,  and  A  yield  the  curves  of  figure  (5.  8). 
In  this  plot,  the  normalized  effective  SNR  (6  /SNR)  is  plotted  versus 
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the  %  of  power  in  the  second  signal  for  a  range  of  values  for  A.  The 
curves  show  that,  for  a  given  element  spacing-delay  factor  t,  the 
system  performance  degrades  as  the  multipath  delay  increases,  and 
degrades  as  a  greater  percentage  of  the  available  signal  power  is  trans¬ 
mitted  via  the  second  (later)  path.  These  results  show  the  same  trends 
as  for  the  single  element  system  of  Case  I.  The  degradation  is  again 
the  result  of  loss  of  signal  energy  concentration,resulting  in  lowered 
input  SNR  due  to  multipath  spreading  of  the  random  signal.  The  multi- 
path  structure  thus  causes  the  same  effects  in  the  array  detection  case 
as  in  the  single  element  case. 

ROC  curves,  corresponding  to  the  curves  for  the  index,  are 
plotted  in  figure  5-9  for  fixed  multipath  power  splitting  (A  =  70%,  B  =  30%) 
and  varying  A  and  in  5-10  for  fixed  A  and  varying  power  splitting 
(,5sAsl.O  ,  0  £  B  £  .  5')  In  both  figures  limiting  "worst"  and  "best" 

cases  are  plotted. 

In  the  above  analysis,  the  effects  of  changing  multipath  structure 
were  investigated  for  fixed  element  spacing.  Since  the  multipath  structure 
is  normally  fixed  by  the  physics  of  the  channel,  mechanism  interacting 
with  the  signal,  basically  uncontrollable,  any  improvement  in  system 
performance  due  to  design  must  come  by  variation  of  the  array  geometry, 
steering,  or  frequency  filtering  /  data  processing. 

The  performance  index  of  (5.  22)  is  plotted  in  figure  5-11  for 
fixed  multipath  structure  (A  =  1.  0)  and  various  values  of  relative  element 
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spacings.  These  curves  show  the  pronounced  effect  of  element  spacing 
on  performance.  It  can  be  seen  from  5-9  that  system  performance  is 
enhanced  considerably  by  increasing  the  separation  of  the  elements. 

The  interpretation  of  the  variation  of  performance  with  element  spacing 
is  the  same  as  for  Case  Study  II,  and  illustrated  in  Figs.  5-6,  5-7.  The 
increased  spacing  produces  a  narrower  directional  spatial  beam  pattern 
with  better  resolution  and  the  effects  of  the  new  combining  and  frequency 
filtering  stages  perform  optimal  noise  and  sidelobe  rejection,  resulting 
in  improved  performance.  Thus  array  detectability  resolution  is 
increased  by  separation  of  the  array  elements.  The  ROC  curves  for  a 
fixed  multipath  structure  are  plotted  in  figure  5-12,  showing  the  effects 
of  array  element  separation  upon  system  error  probabilities.  Thus, 
varying  element  geometry  in  the  form  of  separation  gives  one  a  method 
of  combating  degradation  due  to  multipath. 

A  comparison  of  the  results  of  Case  Study  I  and  Case  Study  III 
are  shown  in  figures  5-13,  5-14,  5-15,  5-16.  The  figures  show  a  com¬ 
parison  of  the  performance  index  of  a  single  element  system  and  a  two 
element  array  for  a  range  of  multipath  situations.  It  will  be  noted  that 
in  each  figure,  regardless  of  A,  the  curve  of  the  system  index  d  for  the 
single  element  lies  below  the  curve  for  the  array  index,  corresponding  to 
an  element  spacing  of  t  =  Z.O.  For  t<2.0,  the  single  element  perform¬ 
ance  index  is  superior  to  the  array  system.  Since  t  was  plotted  relative 
to  the  random  signal  bandwidth  in  the  figures,  it  is  related  to  Wc  by 
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where  m  is  a  positive  number  dependent  upon  element  spacing  and 
geometry.  It  follows  from  examination  of  the  curves  of  figures  5-13, 

5-14,  5-15,  5-16  that  the  array  element  spacing  must  be  greater  than 
2.  0  times  the  reciprocal  random  signal  process  bandwidth  in  order  for 
a  two  element  system  to  perform  more  optimally  than  a  single  element 
system.  For  white  noise  increasing  the  spacing  beyond  2.  0/u)^  produces 
even  greater  improvement  in  performance,  and  the  array  system  is 
clearly  preferable.  If  t<2.0,  the  single  element  receiving  system 
is  preferable.  Thus,  the  conclusion  is  drawn  that  the  key  in  choosing 
array  geometry  (element  spacing)  for  a  multipath  signal  environment  is 
that  the  element  spacing  be  significantly  larger  than  the  reciprocal  signal 
process  bandwidth.  An  improved  array  system  detection  performance 
capability  is  realized  by  chossing  an  element  spacing  with  respect  to  the 
power  spectrum  bandwidth  of  the  incident  signal,  not  with  respect  to 
the  multipath  structure  in  this  case.  The  array  geometry  should  be 
chosen  to  give  "best"  performance  possible  (as  widely  spaced  elements 
as  technically  feasible)  under  the  assumption  that  there  is  no  multipath 
structure  in  the  random  signal.  If  t>2,  0/u>c,  this  insures  that  the  array 
will  give  superior  performance  over  the  single  element.  If  a  multipath 
structure  does  actually  exist,  the  performance  will  degrade  simultaneously 
compatible  amounts  in  both  systems,  as  is  shown  by  comparing  figures 
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5-13  through  5-16  If  the  multipath  structure  is  fixed,  the  major  method 
available  to  the  designer  to  improve  a  system's  performance  is  to 
increase  the  element  spacing;  once  the  maximum  spacing  (dictated 
by  practical  size  and  space  requirements)  is  reached,  and  if  t>2,  0/u)c, 
performance  is  optimal.  This  dependence  results  from  the  separation 
of  array  gain  function  (noise  field  and  steering  effects)  and  the  sign?.l-to- 
noise -ratio  (random  signal  field,  multipath  effects)  in  the  performance 
index.  The  use  of  multiple  element  detection  systems  with  proper 
element  spacing  can  be  used  to  overcome  degradation  in  a  single  element 
system  caused  by  the  multipath.  The  above  assumes  that  desired  spacing 
can  be  achieved  such  that  the  receiving  elements  are  located  in  a  manner 
which  renders  cross-link  noise  correlation  zero. 

In  general,  the  preceding  case  study  has  shown  that 

1<  system  performance  degrades  as  the  multipath  delay 
increases 

2.  system  performance  degrades  as  more  signal  energy 
is  transmitted  via  later  paths 

3.  system  performance  improves  as  the  element  separation 
is  increased 

4.  two  element  system  performance  is  optimal  over  a 
single  element  system  for  array  spacing  greater  than  2,  0 
times  the  reciprocal  bandwidth  of  the  random  signal  in 
any  given  multipath  environment,  assuming  approximately 


a  single  arrival  direction. 


The  above  case  study  establishes  the  utility  of  array /diversity 
receiving  systems  and  illustrates  a  method  of  evaluating  attainable  perform¬ 
ance  capability  over  single  receiving  element  systems. 

5-5.  Case  IV  -  Multi-Element  Array,  Optimum  Geometry,  No  Multipath 

The  above  studies  have  shown  that  the  ’'best"  performance  for  a 
multi-element  array  is  realized  for  an  array  with  widely  spaced  elements 
imbedded  in  a  channel  with  no  multipath.  In  such  cases,  the  performance 
characteristics  are  optimum.  ROC  curves  for  a  single  element,  a  two 
element,  and  a  three  element  receiver  in  such  an  optimum  environment 
are  shown  in  figure  5-17.  In  this  plot,  the  performance  of  multi-element 
detection  systems  is  illustrated,  Curves  for  systems  N>3  lie  above  the 
N  =  3  curve  and  have  superior  performance  characteristics.  Thus 
detection  capability  can  be  enhanced  by  properly  employing  optimized 
multi-element  arrays.  In  general,  properly  designed  space-time  multi¬ 
element  array  detection  systems  offer  more  flexibility  and  superior 
system  performance  characteristics  than  do  systems  employing  a  single 
element. 

5-6.  Summary 

The  series  of  case  studies  has  shown  the  effects  upon  system 
performance  (error  probabilities  and  detectability)  in  terms  of 
1)  multipath  structure  in  the  incident  random  signal(multipath  delay  and 
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multipath  power  division),  2)  array  geometry  and  steering  (element 
spacing),  and  3)  space-time  (spatio-temporal)  correlation  properties 
of  the  incident  random  processes.  In  general,  the  following  conclusions 
can  be  established: 

1.  Pronounced  multipath  structure  in  the  random 
signal  degrades  system  performance  due  to  spreading  of  the 
signal  energy  over  time.  Increasing  delay  and  power  split 
degrade  performance,  with  the  delay  factor  being  the  major 
factor  in  degrading  performance. 

2.  Multi-element  systems  suffer  comparable  degra¬ 
dation  in  performance  resulting  from  multipath,  as  do  single 
element  systems. 

3.  Multi-element  systems  offer  better  performance 
characteristics  than  single  element  systems,  if  the  spacing 
of  the  array  elements  is  greater  thanm  times  the  reciprocal 
bandwidth  of  the  'infcident  random  signal-,  where  m  is  a ‘constant 
dependent  upon  the  number  of  elements,  the  signal  and  noise 
fields,  aftd  the  array  geometry.  For  N  =  2,  an  interferometer, 
m  ~  2.  0. 

4.  The  maximum  performance  gain  possible  for  the  cases 
considered  is  realized  for  large  element  spacing. 

5.  The  array  system  requirements  and  performance 
characteristics  are  highly  dependent  upon  the  spatial  structure 
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of  the  array  and  the  space-time  structure  of  the  incident 
signal  and  interference  fields, 

6.  The  array  space-time  analysis  technique  is  an 
extremely  powerful  tool  for  analysis,  design,  implementation, 
and  evaluation  of  modern  detection  systems. 


CHAPTER  VI 


CONCLUDING  REMARKS  AND  FUTURE  RESEARCH 

It  has  been  the  purpose  of  this  research  effort  to  obtain  definitive 
solutions  to  the  problems  associated  with  the  analysis  and  design  of  optimal 
space(array)diversity  receiving  systems.  In  particular,  it  is  of  consider¬ 
able  technological  interest  to  establish  the  optimum  system  requirements 
and  system  configurations,  and  to  ascertain  their  related  performance 
characteristics  for  a  wide  range  of  physical  information  channel,  cross- 
section  and  noise  field  environments.  The  results  and  conclusions 
obtained  in  this  effort  offer  considerable  design  improvement  capability 
in  the  areas  of  sonar /radar /communications/ seismology,  as  well  as 
offering  a  unified,  more  concise  formulation  for  many  general  reception/ 
detection/estimation  problems. 

This  effort,  along  with  other  studies  in  spatio-temporal  reception 
(Refs.  1,  2, 14,  20-24,46)  has  shown  that  the  dichotomy  of  the  detection 
problem  into  tv/o  separate  areas  (antenna  and  sensor  design,  and  statis¬ 
tical  communication  theory)  can  result  in  design  and  implementation  of 
receiver  systems  that  are  not  truly  optimum,  in  the  sense  of  maximum 
exploitation  of  the  available  information  and  physics  of  the  environment. 
Space-time  considerations  attempt  to  combine  the  spatial  properties  and 
the  temporal  properties  of  the  incident  processes  physics,  the  array, 
its  geometry,  and  the  associated  data  processor  in  order  to  obtain 
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systems  that  are  optimum  in  both  their  space  and  their  time  data  process¬ 
ing  characteristics.  It  has  been  shown  that  such  systems,  when  properly 
designed,  offer  performance  improvement  over  sub-optimum  systems 
implemented  with  limited  "classically  optimal"  procedures.  This 
improved  performance  is  not  obtained  without  cost,  however,  since  the 
analysis  and  design  of  space-time  systems  is  considerably  more  difficult 
than  conventional  methods.  Thus,performance  enhancement  is  possible  by 
the  use  of  spatio-temporal  techniques,  enhancement  that  is  paid  for  by 
larger,  more  complex,  yet  quite  physically  tractable,  systems. 

Finally,  the  effort  has  shown  that  in  the  final  analysis,  in  design 
and  implementation  of  receiving  systems,  the  engineer  must  deal  with 
the  real,  basic  physics  of  wave  mechanics  and  receiving  element- wave 
front  interaction.  The  mathematical  elegance,  abstractness,  and  com¬ 
plexity  that  arises  in  the  analysis,  is  only  of  value  if  it  reflects  and 
models  the  actual  physical  mechanisms  that  exist  in  the  physical  environ¬ 
ment  of  the  systems.  The  goal  in  engineering,  and  in  particular,  in 
space-time  reception,  is  to  express  the  system  requirements  and  resulting 
performance,  in  terms  of  the  fundamental  variables  of  physics  (space, 
time  and  energy).  Solutions  not  in  these  terms,  or  not  capable  of  being 
reduced  to  them,  are  of  little  practical  engineering  value.  It  has  been 
the  purpose  of  this  research  to  express  the  results  in  terms  of  these 


variables. 


This  research  effort  has  quantitatively  established  the  optimum 
system  requirements  and  configurations,  and  quantitatively  established 
the  effects  upon  system  performance  of 

1.  multipath  structure  in  the  random  signal,  both 
variations  in  the  path  delay  and  power  division  between  thevarious  paths. 

2.  the  effects  of  array  element  separation  (array 
geometry)  relative  to  the  struitu^e  of.the  incident  processes 

3.  the  effects  of  certain  spatio-temporal  correlation 
properties  of  the  incident  processes. 

In  particular,  it  was  found  that 

1.  increasing  multipath  delay  degrades  performance,  as 
dees  spreading  of  the  available  signal  power  over  a  number  of 
transmission  paths 

2.  in  array  detection,  the  key  relationship  is  the  signal 
bandwidth  to  element  spacing.  If  the  spacing  is  such  that  it  is 
approximately  2.  0  times  the  reciprocal  bandwidth  of  the  random 
signal,  the  array  will  perform  in  a  superior  manner  when  com¬ 
pared  to  a  single  element  system.  Decreasing  the  spacing  below, 
this  value  yields  the.  single  receiving  element  system  preferable. 

3.  the  optimum  system  configuration  in  a  two  element 
array  is  such  that  the  elements  are  as  widely  separated  as 
possible. 


Finally,  some  criticism  of  the  preceding  analysis  is  necessary; 
these  criticisms  give  rise  to  future  necessary  studies  and  further 
research  efforts. 

The  most  obvious  need  is  for  the  analysis  to  be  extended  to  multiple 
element  (N>  2)  arrays.  The  analysis  methods  presented  here  are  valid 
for  these  larger  arrays.  This  study  (N  =  2)  gives  insight  into  the  effects 
and  shows  the  general  trends  of  the  effects  of  element  spacing  and  multi- 
path  upon  system  performance.  Analysis  for  (N>  2)  will  give  more 
definitive  answers  as  to  the  effect  of  element  spacing  and  geometry  in 
multi-element  systems. 

Additional  study  is  needed  for  space-time  random  fields  other 
than  those  considered  in  the  case  studies.  In  particular,  the  effects  of 
the  magnitude  of  cross  link  correlation  of  the  noise  and  signal  fields 
need  to  be  established.  The  analysis  presented  above  considered  only 
white  noise  and  Markovian  signals;  other  distributed  noise  field  spectra 
should  be  investigated. 

The  multipath  arrival  angle  in  the  cases  was  assumed  essentialy  the 
same  for  all  paths.  In  practice,  the  angles  are  usually  different;  this 

46 

analysis  needs  to  be  investigated  further.  A  recent  effort  by  Schweppe 
has  initiated  a  study  of  this  area. 

18 

The  signals  were  considered  purely  random.  Middleton  has 
shown  that  a  strong  specular  component  is  present  quite  often  in  the  cases 
of  random  fading,  scattering,  and  reverberation  of  deterministic 


signals.  This  result  should  be  incorporated  into  the  system,  design 
equations. 

Finally,  the  analysis  was  carried  out  by  assuming  that  the  coup¬ 
ling  was  memoryless.  The  analysis  should  include  the  specific  effects 
of  the  interaction  of  the  receiving  elements  and  the  space-time  wave 
equations.  Such  multivariable  analysis^’  ^  would  strengthen  the  argu¬ 
ments  for  consideration  of  the  unified  space-time  reception  problem. 

In  effect,  this  research  is  only  the  beginning  of  a  much  broader, 
more  general  research  effort  that  has  been  critically  needed  in  studying 
the  reception  problem.  There  is  a  definite  need  to  reorganize  and  unify 
the  disciplines  of  statistical  communication  theory  and  antenna  array 
design  theory  in  order  to  achieve  a  unified  approach  to  the  solution  of 
the  reception  problem.  This  research  effort  has  investigated  a  small 
area  of  the  total  problem.  It  has  shown  that  space-time  analysis  is 
definitely  valuable  in  system  design  and  has  given  definitive,  quantitative 
answers  to  the  case  study  problems  of  optimum  array  detection  of  noisy 
random  multipath  signals.  It  has  established  the  effects  of  multipath 
(delay  and  power  division)^  of  element  geometry  (spacing^  and  of  spatio- 
temporal  correlation  properties  of  the  incident  random  fields  on  performance 

The  over -all  result  of  this  research  effort  is  the  establishing  of 
general  synthesis  guidelines  that  enable  one  to  attain,  design,  implement^ 
analyze^and  enhance  performance  of  optimum  space-time  detection/ 
estimation/classification/  and  discrimination  systems. 
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characteristics  of  maximum-likelihood  space-diversity  receiving  systems.  From  the  maximum- 
likelihood  principle,  and  the  multidimensional  Karhunen-Loeve  expansion,  the  continuous  decision 
equations  for  space- time  detection  of  noisy,  random  multipath  signals  are  developed.  Interpretation 
of  these  equations  yields  the  optimum  space-time  receiving  system  configuration  requirements. 

The  performance  characteristics  of  the  resultant  systems  are  then  analyzed  to  determine  the  effects 
of  channel  multipath  structure  (multipath  delay  and  power  division  among  the  paths),  space-time 
correlation  properties  of  the  Incident  processes,  and  the  specific  temporal  correlation  introduced  by 
the  array  geometry.  It  is  shown  by  a  series  of  case  studies,  that  for  both  the  single  element  coupling 
as  well  as  array  multiple  element  coupling,  that  increasing  the  multipath  delay  factor  results  in 
decreased  system  performance  capability  for  fixed  power  of  the  signal  and  noise  processes.  Similarly, 
the  performance  capacity  is  degraded  as  the  available  signal  power  tends  to  distribute  more  equally 
over  the  multiple  transmission  paths.  These  effects  are  attributed  to  the  loss  of  effective  signal 
energy  concentration,  resulting  In  a  lowei  effective  detectabilUy  signal-to-noise  ratio.  An  investiga¬ 
tion  of  these  effects  upon  system  performance,  due  to  array  geometry  (receiving  element  spacing) 
shows  that  performance  is  enhanced  by  increasing  the  separation  distance  of  the  elements  for  a  given 
multipath  situation.  This  phenomena  is  the  result  of  a  more  compatible  array  beam  pattern  that  comes 
■about  from  Increased  spacing.  In  effect,  with  increased  spacing,  the  main  lobe  of  the  pattern  is 
narrowed,  while  the  side  lobes  are  optimally  suppressed  by  the  optimal  cross-coupling  of  elements. 
Finally  the  analysis  illustrates  how  optimum  space-time  receiving  systems  capability  results  from  a 
joint  consideration  of  the  coupling  and  temporal  processing,  rather  than  a  dlchotomization  of  the 
problem  into  a  spatial  (antenna)  problem  and  a  (temporal)  data  processing  problem.  As  such,  the 
analysis  yields  a  more  concise,  broader,  interpretation  of  system  design  requirements  and  associated 
performance  characteristics. 
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